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Title: Method of Calculating Shading Correction Coefficients of Imaging Systems 
from Non-Uniform and Unknown Calibration Standards 

[01] This application claims the benefit of U.S. Provisional Patent application 

No. 60/236,945, which is incorporated herein by this reference. 

Field of the Invention 

[02] The present invention relates to a process for analyzing shading errors in 

an imaging system and determining shading correction coefficients to compensate for 
non-uniformity in illumination, optics and imaging sensors and electronics. 

Background of the Invention 

[03] Images acquired by an imaging system are not perfect. One common 

type of distortion is shading error (systematic intensity variation across acquired 
images), which is mainly caused by: 

i. Photo-electronic conversion (sensors) inhomogenity; 

ii. Illumination inhomogenity; 

iii. Optical inhomogenity; 

iv. Electronic inhomogenity. 

[04] If not corrected, shading error will distort quantification performed from 

image intensities. It will also distort the visual appearance of the acquired images. 
[05] In the prior art, shading errors and shading correction coefficients of 

imaging systems are calculated from the image of a calibration standard (or reference) 
with uniform reflectance, transmittance, or emission across the entire field of view. 
Typically, the materials in a calibration standard are the same or have the same 
characteristics of the materials to be imaged. For example, to correct a fluorescent 
imaging system, a slide covered by a thin layer of fluorescent material with uniform 
thickness and density may be used as a calibration standard. The fluorescent 
substance in the calibration slide should have the same excitation/emission properties 
as those to be imaged and/or analyzed subsequently. 



[06] 



Typically, an imaging system is used to image and/or to analyze objects in 



a continuous fashion. That is, every point in the field of view is of interest. In this case, 
pixel-by-pixel shading corrections should be performed. Figure 1 illustrates the 
mapping between object plane and image sensor plane. The rectangular area 200 in 
the object plane of an imaging system 250 is projected to the image sensor 202 of the 
imaging system 250. The area 200 is also referred to as the field of view of the system 
250. In the example shown in Figure 1 , the image sensor 202 consists of 48 sensor 
elements 206 (or pixels) in an 8 x 6 array, denoted as pxl(x, y), x=1,..,8 and y=1,...,6. 
Correspondingly, there are 8 x 6 divisions 205, denoted as Rp(x, y), x=1,...,8 and 
y=1,...,6, in the field of view 200 in the object plane. The region Rp(x, y) is projected 
and focused to sensor element pxl(x, y) through the lens 201 of the imaging system 
250. In this example, there will be 8 x 6 shading correction coefficients, each of which is 
an integrated (average) coefficient of the corresponding region. 
[07] If the image sensor of an imaging system has N P by Mp elements, and 

/C P (x, y) is the image pixel intensity of pixel pxl(x, y) of an uniform calibration standard, 
the pixel-by-pixel shading error SE P (x, y) is defined in EQ. 1. Typically, shading errors 
are normalized so that the average of all errors equals to one. 



levels of the imaging system have a linear relationship with the input signal. For 
example, when imaging an illuminant object, the input signal is the self-emitted light 
intensity from the object. In the linear case, image intensity /C P (x, y) is the raw pixel 
gray-level of the pixel pxl(x, y) minus a stored background value of the corresponding 
pixel. In the non-linear case, the non-linear relationship first is modeled, and then this 
non-linear model is used in calculating the imaging intensity /Cp(x, y) from the raw gray- 
levels so that the image intensity IC P {x, y) is linear with the input signal. If the image 
intensity is not linear to the input signal or the background values are not properly 
removed, the calculated shading error will not be accurate. In such a system, a skilled 



SE p (x,y) = 




Forx=1,...,N P and y=1,...,M P EQ. 1 
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Typically, an imaging system is a linear device. That is, the raw pixel gray 



person will be capable of adjusting the image intensity values so that they have a linear 
relationship with the input signal. 

[09] Figure 2 generalizes the pixel-by-pixel shading error concept. Instead of 

establishing shading correction coefficients of a region corresponding to each single 
pixel for the entire field of view, shading correction coefficients can be established for 
regions corresponding to more than one pixel. The pixel-by-pixel example of Figure 1 is 
a special case of this generalized model. 

[10] In the example shown in Figure 2, there are 3 x 2 regions 215, denoted as 

R(x, y): x=1 ,2,3 and y=1 ,2. Each of these regions projects to its corresponding pixel set 
216, denoted as P(x, y): x=1 ,2,3 and y=1 ,2, in the imaging sensor 202. Each pixel set 
in this example contains four pixels. For example, P(1 , 1), corresponding to R(1 , 1), 
consists of: {pxl(1, 1), pxl(2, 1), pxl(1, 2), pxl(2, 2)}. The generalized shading correction 
has 3x2 coefficients, and each coefficient corresponds to an area of 4 pixels. Of 
course, the regions do not have to be in a regular 2-dimensional matrix form. They can 
be in any form. For the convenience of discussion, the regular 2-dimensional matrix 
form is used in this document. The calibration standard may be 1 -dimensional or 3- 
dimensional or may have a non-rectangular 2-dimensional form. 
[11] Given a calibration standard with uniform regions R(x, y): x=1,...,N and 

y=1,...,M, the shading error SE(x, y) for the region R(x, y) (or pixel set P(x, y)) is defined 
in EQ. 2. 

SE(x 9 y)= IC N (x > M y) Forx=1,...,Nandy=1,...,M EQ, 2 

where: /C(x, y) is the image intensity of region R(x, y). It is calculated by 

averaging the image pixel intensities of all the pixels corresponding to 
region R(x, y). 

[12] Once the shading error is defined, the shading correction coefficient SC(x, 

y) for the region R(x, y) is defined in EQ. 3. 

SC(x,y) = — - — Forx=1,...,N andy=1,...,M EQ, 3 

SE(x 9 y) 

[13] The shading correction coefficients SC(x, y) will then be stored and 

applied to subsequent images to compensate for the system's shading error. 
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[14] The key requirement of the prior art is uniformity. However, uniform 

calibration standards are often difficult or impossible to prepare, and, therefore, the 
resulting estimated shading correction coefficients are often inaccurate. 
[15] Another method for calculating the shading correction coefficients is using 

calibration standards with non-uniform but known local reflectance, transmittance or 
emission values. Again, it is often very difficult to accurately determine the required 
reflectance, transmittance or emission values at all the required locations, especially in 
microscopic imaging applications. 

Summary of the Invention 

[16] The present invention provides a method for calculating the shading 

correction coefficients for an imaging system from non-uniform and unknown calibration 
standards. 

[17] A calibration standard which corresponds to the item or material to be 

imaged is first selected. The calibration standard may have one, two or three 
dimensions. Potentially, the calibration standard may have more dimensions. 
Preferably, the calibration standard covers all or most of the field of view of the imaging 
system. The calibration standard is divided into regions and imaged to calculate an 
initial image intensity. The calibration standard is re-arranged and imaged again to 
determine a set of re-arranged image intensities. The process of re-arranging and 
image is repeated once for each dimension of the calibration standard. 
[18] The initial image intensity and re-arranged image intensities are combined 

mathematically to form a set of equations in which a shading error for each region may 
be calculated based on the shading error of another region or regions. By selecting a 
value for the shading error of one region, an intermediate shading error of each region 
may be calculated. The shading error may then be calculated as the normalized value 
of the intermediate shading error. The shading correction coefficient for each region 
may be calculated by inverting the shading error for that region. 
[19] The shading correction coefficients may then be used to correct the image 

intensity values for an actual object imaged using the imaging system. 



[20] The calibration standard need not have uniform regions, as long as each 

region has the same characteristics as the object to be imaged. Differences between 
the regions are abstracted out of the determination of the shading correction coefficient 
by the present invention. 

[21] In one aspect, the present invention provides a method for calculating 

shading correction coefficients for an imaging system, comprising the steps of defining a 
set of calibration regions on a calibration standard; defining a first part and a second 
part of said calibration standard, wherein each of said first and second parts contains at 
least one calibration region; taking an image of said calibration standard and calculating 
an initial image intensity for each of said calibration regions; re-arranging said first and 
second parts to form a re-arrangement of said calibration standard; taking an image of 
said re-arrangement and calculating a re-arranged image intensity for each of said 
calibration regions; calculating a shading error for each of said calibration regions; and 
calculating a shading correction coefficient for each of said calibration regions. 
[22] In another aspect, the present invention provides a method for calculating 

shading correction coefficients for an imaging system, comprising the steps of: defining 
a set of calibration regions on a calibration standard; defining a first part, a second part, 
a third part and a fourth part of said calibration standard, wherein each of said first, 
second, third and fourth parts contains at least one calibration region; taking an image 
of said calibration standard and calculating an initial image intensity for each of said 
calibration regions; re-arranging said first, second, third and fourth parts in a first 
direction to form a first re-arrangement of said calibration standard; taking an image of 
said first re-arrangement and calculating a first re-arranged image intensity for each of 
said calibration regions; re-arranging said first, second, third and fourth parts in a 
second direction to form a second re-arrangement of said calibration standard; taking 
an image of said second re-arrangement and calculating a second re-arranged image 
intensity for each of said calibration regions; calculating a shading error for each of said 
calibration regions; and calculating a shading correction coefficient for each of said 
calibration regions. 
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Brief Description of the Drawings 

[23] Several exemplary embodiments of the present invention will now be 

described, by way of example only, with reference to the drawings, in which: 

Figure 1 is an illustration of pixel mapping of an imaging system; 

Figure 2 is an illustration of generalized region mapping of an imaging 

system; 

Figure 3 is an illustration of a calibration standard with a top-left 

configuration; 

Figure 4 shows the dimensions of the calibration standard in Figure 3; 

Figure 5 is an illustration of a horizontally shifted-and-rotated arrangement 
of the calibration standard of Figure 3; 

Figure 6 is an illustration of a vertically shifted-and-rotated arrangement of 
the calibration standard of Figure 3; 

Figure 7 is a flowchart illustrating procedures for establishing shading 
correction coefficients of an imaging system using the calibration standard of Figure 4; 

Figure 8 is an illustration of a calibration standard with a bottom-right 

configuration; 

Figure 9 is an illustration of a horizontally shifted-and-rotated arrangement 
of the calibration standard of Figure 8; 

Figure 10 is an illustration of a vertically shifted-and-rotated arrangement 
of the calibration standard of Figure 8; 

Figure 11 is an illustration of a calibration standard with parallelogram 

form; 

Figure 12 is an illustration of a calibration standard with a top-left 
configuration and with two calibration regions in a corner part; 

Figure 13 is an illustration of calibration standard with a top-left 
configuration, but which cannot be disassembled; 

Figure 14 illustrates a method of performing two shifts to imitate a 
horizontally shifted-and-rotated arrangement of the calibration standard of Figure 13; 

Figure 15 illustrates a method of performing two shifts to imitate a 
vertically shifted-and-rotated arrangement of the calibration standard of Figure 13; 



Figure 16 is an illustration of a one-dimensional calibration standard with a 
left configuration; 

Figure 17 shows the dimensions of the calibration standard of Figure 16; 

Figure 18 is an illustration of a shifted-and-rotated arrangement of the 
calibration standard of Figure 16; and 

Figure 19 is an illustration of a three-dimensional calibration standard with 
a top-left-front configuration. 

Detailed Description of the Invention 

[24] The present invention provides a method for calculating shading 

correction coefficients for an imaging system, using a calibration standard with N x M 
regions. Reference is first made to Figure 3, which illustrates a first calibration standard 
10 according to the present invention. The described calibration standard 10 consists of 
N x M (N columns and M rows) calibration regions 5, which are referred to as R(x, y). 
Calibration standard 10 is divided into four parts. The variables x and y are the array 
index, x takes values form 1 to N, denoted as x = 1,...,N, and y takes values from 1 to 
M, denoted as y = 1 ,M. The top-left corner 4 consists of one calibration region 5 at 
its center; the top-row part 2 consists of N-1 calibration regions 5; the left-column part 3 
consists of M-1 calibration regions 5, and bottom-right part 1 consists of (N-1) x (M-1) 
calibration regions 5. In this exemplary embodiment, N > 1 and M > 1 . This is referred 
to as the 2D (2-dimensional) case. In the 2D case, both N and M have to be greater 
than one. If either N or M is equal to one, it becomes a one dimensional case, which is 
described later later. 

[25] The location of part 4 (i.e., the part with one calibration region R(1 ,1)) of 

calibration standard 10 is used to name the configuration of the calibration standard. 
The calibration standard shown in Figure 3 has a top-left configuration, and the 
arrangement of the four parts shown in Figure 3 is the initial setup of the top-left 
configuration. 

[26] Calibration standard 10 is imaged by an imaging system and a shading 

correction coefficient for each region 5 is calculated according to the present invention. 
The shading correction coefficients for regions 5 will depend on the object being imaged 
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and the particular characteristic being imaged. For example, a fluorescent object may 
be imaged to determine the amount of excited light from the object Subsequently, the 
same imaging system may be used to image an illuminant object to determine the 
amount of self-emitted light. In these two different cases, the shading correction 
coefficients may be different, even for comparable regions of the imaging field. In the 
first case, the shading error is due to the inhomogenity of the lens, the excitation light 
source, sensor and other electronic sources. In the second case, the shading error is 
due to the inhomogenity of the lens, sensor and other electronic sources (there is no 
excitation light source). 

[27] It is preferable that the material in each calibration region 5 is the same or 

has characteristics similar to the object which will ultimately be imaged. For example, 
the same fluorescent material is preferably deposited on the calibration regions 5 if the 
object to be imaged is fluorescent. However, unlike the prior art, it is not necessary that 
the densities or concentrations of the deposited material in regions 5 of calibration 
standard 10 are the same and the differences between regions do not need to be 
known. 

[28] The physical shape of the calibration regions 5 can be in any form. The 

square shape shown in the Figures is for illustration purpose only. Also, the gaps 
between calibration regions 5 are for the purpose of illustration only. Their physical 
presence is not required. Also, the central placement of calibration regions 5 in parts 1, 
2, 3 and 4 are for the illustration purpose only. Their placement could be off the center 
as long as the region intervals remain uniform. 

[29] The signal intensity of region R(x, y) is denoted as VC(x, y). In the case of 

illuminant calibration standards, the signal intensity is the emitted light intensity. In the 
case of fluorescent calibration standards, the signal intensity is the excited light intensity 
under uniform illumination since the illumination inhomogenity is already incorporated 
into the shading error. The signal intensity VC{x, y) of a particular region R(x, y) should 
remain constant during the calibration process described below. 
[30] The dimensions of the calibration standard 10 are shown in Figure 4, in 

which "a" is the horizontal interval of regions 5 and "b" is the vertical interval of regions 
5. Preferably, all the calibration regions 5 should be within the field of view of the 
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imaging system. Also, preferably, the spatial dimension of the calibration standard 10 
should be equal to that of the field of view of the imaging system to be calibrated. That 
is: 

Wfov = a x N and Hfov ~ bxM 

where: Wfov and Hfov are the width and the height of the field of view of the 
imaging system to be calibrated. 
[31] The selections of a, b, N, M and the shape of the calibration regions 5 are 

determined by requirements of a particular application of the present invention. Where 
a large variation in shading errors is expected, a and b may be selected to be small and 
consequently N and M (assuming a fixed field of view) will be large. That is, the regions 
5 are spaced more closely together. 

[32] When the calibration standard shown in Figure 3 is imaged in an imaging 

system, the imaging intensity/C°(>,.y) of the region R(x, y) is defined by the following 

equation: 

IC°(x 9 y) = GxSE(x,y)xVC(x 9 y) EQ. 4 

forx=1,...,N and j/=1,...,M 

where: G is a constant representing the imaging system gain. 

[33] Taking a natural log transform on both side of EQ. 4, we obtained: 

IC°(x 9 y) = G L +SE L (x 9 y) + VC L (x 9 y) EQ. 5 

for x=1,...,N and y=1,...,M 

where: IC° L (x 9 y) = In (IC°(x 9 y) ) 

SE L (x, y) = In (SE(x, y)) 
VC L (x, y) = In (VC(x, y)) 
G L = In (G) 

In (d) is the natural log of number d. 
[34] Figure 5 shows a horizontally shifted-and-rotated arrangement 1 1 of the 

top-left configuration of calibration standard 10 (Figures 3 and 4). The horizontal shift- 
and-rotation is done by shifting the entire calibration standard 10 to the left by one 
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column and then moving parts 3 and 4 to the right of parts 1 and 2. The re-arranged 
calibration regions 5 are referred to as R H (x, y), and the relationship between the 
contents (calibration materials) of R H (x, y) and R(x, y) is defined as: 

Content(R* (jc,y)) = Content(R(M<wf (jc + lN\y)) EQ. 6 

for x=1,...,N and y=1,...,M 

[35] The function Mod (q t d) in EQ. 6 and subsequent equations is defined as: 

(q + mxd-Y)%d + l if q and d are integers and q . < landd>l 
Mod(q,d) ~) (q~ 1)% d + l if q and d are integers and q > 1 and d > 1 

L Invalid Otherwise 

where: % is a remainder operator (e.g., 1%3 = 10%3 =1). 

Given integer q and d, m is the smallest positive integer such that q + mx 

d> 1. 

[36] Figure 6 shows a vertically shifted-and-rotated arrangement 12 of the top- 

left configuration of calibration standard 10 (Figures 3 and 4). The vertical shift-and- 
rotation is done by shifting the entire calibration standard 10 up by one row and then 
moving parts 2 and 4 to the bottom of parts 1 and 3. The re-arranged calibration 
regions 5 are referred to as R v (x, y), and the relationship between the contents of R v (x, 
y) and R(x, y) is defined as: 

Content(R v (x 9 y)) = Content(R(x,M^( y + 1,M))) EQ. 7 

for x=1,...,N and y=1,...,M 

[37] According to the present invention, the re-arranged calibration regions 5 of 

the shifted-and-rotated arrangements 11 and 12 are physically aligned to the calibration 
regions 5 in the initial setup of calibration standard 10. That is: 

P(x, y) = P H (x, y) = P v (x, y) RQ. 1 

for x=1 ,...,N and y=1,...,M 

where: P(x, y), P H (x, y), and P v (x> y) are the mapped pixel sets of calibration 
regions R(x, y) (Figure 3), R H (x, y) (Figure 4) and R v (x> y) (Figure 5) 
respectively. 
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[38] Let IC?(x,y) be the natural log transformed image intensity of the region 

R H (x, y) for the horizontally shifted-and-rotated arrangement 11 shown in Figure 5. 
From RQ. 1 , we know that the region R(x, y) and the region R H (x, y) map to the same 
pixels. Therefore, they share the same shading error SE(x, y). From EQ. 6, we further 
know that the content of region R H (x, y) is the content of region R(Mod(x+1 , A/), y), i.e., 
the signal intensity of calibration region R H (x, y) is \/C(Moc/(x+1 , A/), y). Therefore, we 
have: 

/Cf (x, y) = G L + SE L (x, y) + VC L (Mod(x + 1, N), y) EQ. 8 

for x=1,... J N and y=1,...,M 

[39] Let IC v L {x,y) be the natural log transformed image intensity of the region 

R v (x, y) for the vertically shifted-and-rotated arrangement 12 shown in Figure 6. Similar 
to the deduction of EQ. 8 from EQ. 7 and RQ. 1, we have: 

IC V L (x, y) = G L + SE L (jc, y) + VC L (x, Mod(y + 1 ? M)) EQ. 9 

for x=1 I ...,N and y=1,...,M 

[40] Multiplying 2 to both sides of EQ. 5, multiplying -1 to both sides of EQ. 8 

and EQ. 9, and then adding the three new equations together, i.e., 2 x (EQ. 5) - (EQ. 8) 
- (EQ. 9), we have: 

/CV(x,;y) = 2xVC L (x,y)-Vq^ EQ. 10 

for x=1,...,N and y=1,...,M 
Where: ICV(x,y) = 2x IC° L (x 9 y)- ICl (x 9 y) - ICl (x, y ) 

[41] Using EQ. 5, we can obtain the following relationship: 

VC L (x,y) = IC° L (x,y)-G L -SE L (x,y) 

VC L (Mod(x + l,A0,y) - IC° L {Mod(x + l,Nly) - G L - SE L (Mod(x + l,N\y) 
VC L ( jc, Mod(y + 1, M)) = IC° L (x, Mod(y + 1, M)) — G L — SE L (x, Mod(y + l,Af)) 

[42] Substituting VC L (x,y) , VC L (Mod(x + l,N),y) and VC L (x,Mod(y + 1,M)) in 

EQ. 10 with the above relationship provides EQ. 11: 
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ICS(x,y) = 2xSE L (x,y)-SE L (Mod(x + l,N),y)-SE L (x,Mod(y + l,M)) EQ. 11 
forx=1,...,N and y=1,. ..,M 



Where : ICS(x, y) = IC H L {x,y) + IC\ (x,y)- IC° L (Mod(x + \,N),y)- IC° (x,Mod(y + 1, M)) 

[43] If we order the two dimensional index (x, y) in EQ. 1 1 by the column first 

method, i.e., (1, 1), (2, 1), (N, 1), (1,2), (2, 2), (N, 2), (1, M), (2, M), (N, 
M), we can represent EQ. 1 1 in a matrix representation, as defined in EQ. 12. 

ICS = A SE|_ EQ. 12 

where: A is a (NxM) x (NxM) matrix, i.e., it has NxM rows and NxM columns. It 
can be represented by a partitioned matrix, as defined below. The 
matrix E is a N x N matrix as defined below, the matrix I is a N x N 
identity matrix and matrix 0 is a N x N zero matrix (all coefficients are 0). 



A = 



E -I 
0 E 
0 0 



0 

-I 

E 



0 0 
0 0 
-I 0 



0 
0 
0 



0 0 
0 0 
0 0 



E -I 0 
0 E -I 
0 0 E 



NxM xNxM 



E = 



2 -1 
0 2 



0 0 
-1 0 



0 0 
0 0 



2 -1 
0 2 



1 = 



NxN 



1 0 

0 1 

0 0 

0 0 



0 0 

0 0 

1 0 
0 1 



NxN 
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ICS = 



7CS(U) 
ICS{%\) 

ICS{N,\) 



ICS(\,M) 
ICS(2,M) 

llCS(N,M)\ x 



SE L : 



(NxM) 



SE L (1,1) 
SE L (2,1) 

SE L (N,l) 



SE L (l,M) 
SE L {2,M) 

SE L (N,M). 



\x(NxM) 



[44] Using standard mathematical techniques, it can be proved that the rank of 

matrix A is NxM » 1 . Therefore the solution of S£t(x, y) is not unique. However, since 
the rank of matrix A is NxM - 1 , SEt(x, y) can be resolved by arbitrarily specifying only 
one of the values (e.g., SE/.(1, 1) is specified). 

[45] Let SEO be a solution of EQ. 11 or EQ. 12, by giving one of the shading 

errors a known value. SEO can be defined as: 

{SEO(x,y,ax,ay,K): x=1 ,...,N, y=1,...,M, 1 < ax < N, 1 < ay < M, DF. 1 

K is a given value and SE L (ax,ay) = K} 

Taking the exponential transform (inverse of the log transform) and normalization (as in 
EQ. 2), the corresponding shading error SE(x, y) is defined as: 

SEO(x,y,ax,ay,K ) 



SE(x,y) = 



a N M 

1 yy 

NxM^ 



EQ. 13 



SEO(i,j,ax,ay,K) 



for x=1,...,N and y=1,...,M 

[46] Using standard mathematical techniques, it can be proved that the 

difference between any given two solutions SEO(x,y,axl,ayl,Kl) and 
SEO(x,y,ax2,ay2,K2) of EQ. 12 is constant, for a given set of measured image 
intensities ( icf(x 9 y) , IC?(x,y) and/C[(x,j;) ). That is: 
SEO(x, y, ax2, ay2, K2) - SEO(x, y 9 axl 9 ayl 9 Kl) = F 
for x=1,... I N and y=1,...,M 
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where: F is a constant. 



[47] 



This give the following equation: 



SEO(x,y,ax2,ay2 f K2) 



e 



SEO(x,y,cucI,ayl,Kl)+F 



-i N M 

NxM^ 



SEO(i,j 9 ax2,ay2,K2) 



i # M 



SEO(iJ,axl,ayl y Kl)+F 



iVxM 
e 



i=\ 7=1 

SEO(x,y,axl,ayl y Kl) 



,SEO(i,j,axl,ayl,Kl) 



[48] That is, the shading error SE(x, y) defined in EQ. 13 is unique, regardless 

of the values of ax, ay and K. 

[49] As an example, select ax=N, ay=M and K=1 as the solution for calculating 

the shading error. This solution is referred to as intermediate shading error and 
denoted as S1(x, y), i.e., Sl(x,y) = SEO(x, y,N,M,l) , whose solution is defined by EQ. 
14Aand EQ. 14B. 

EQ. 14A 



Sl(N,M) = l 
/CS(1,1) 
/C5(2,l) 

/CS(AU) 



ICS(N,M-D + 1 
ICS(l,M) 
ICS(2,M) 



=A1 



57(1,1) 
57(2,1) 

S1(N,V) 



57(1, M) 
57(2, M) 

S1{N-\,M)\ X 



EQ. 14B 



Ix(WxM-l) 



ICS(N-lM) + l_ 

where: Matrix A1 is the top-left (NxM-1) x (NxM-1) sub-matrix of A defined in 



EQ. 12. That is: 



31 lj ~~ 3ij 



for /=1,...,A/-f and y=1,...,M-l 



a1 jj and a,j are the coefficient of matrix A1 and A at column / and row y 
respectively. 
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[50] The constant "1" in the terms "ICS(N, M-1) + 1" and "ICS(N, M-1) + 1" 

results from moving the variable SE L {N, M) to the left side of the equation and SE L (N, 
M) = 1. 

[51] EQ. 14B can be solved by directly calculating the inverse matrix A1" 1 , as 

defined in EQ. 15. It can also be solved by standard iterative methods. 

/C5(l,l) 
/C5(2,l) 



51(1,1) 
5/(2,1) 



5/(AU) 



5/(1, M) 
5/(2, M) 

Sl(N-l,M). 



= A1 



ICS(N,\) 



ICS(N,M -1) + 1 
/C5(1,M) 
/C5(2,M) 

/C5(iV-l,M) + l 



EQ. 15 



lx(ATxM-l) 



[52] The shading error SE(x, y) can then be calculated from the intermediate 

shading error by EQ. 16: 



SE(x 9 y) = 



! N M 

1 yy e si(i 

NxM 



EQ. 16 



j) 



for x=1,...,N and y=1,...,M 

[53] Once shading error SE(x, y) is determined, the shading correction 

coefficients SC(x, y) can be determined from EQ. 3. 

[54] Figure 7 summarizes the shading correction coefficient calculation. A 

calibration standard 10 with a top-left configuration is present to an imaging system to 
be calibrated. Step 101 utilizes the calibration standard 10 in its initial setup as 
illustrated in Figure 3. Step 102 manipulates calibration standard 10 to form a 
horizontally shifted-and-rotated arrangement 1 1 (Figure 5). Step 103 manipulates the 
standard to form the vertically shifted-and-rotated arrangement 12 (Figure 6). In step 
104, images are taken of initial setup of calibration standard 10, shifted-and-rotated 
arrangement 11 and shifted-and-rotated arrangement 12, and the pixel sets P(x, y), 
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P H (x> y) and P v (*. Y) of tne corresponding calibration regions R(x, y), R H (x, y) and R v {x, 
y) are determined using standard image processing algorithms. Step 105 calculates the 
image intensity IC°(x, y), IC H (x, y) and IC v {x, y) from the pixel sets P(x, y), P H (x, y) and 
P v {x, y) and the acquired images. Step 106 calculates the intermediate shading error 
S1{x, y) from EQ. 14A and EQ. 14B. Step 107 calculates the shading error SE(x, y) 
from the intermediate shading error S1(x, y), by EQ. 16. Finally, step 108 calculates the 
shading correction coefficient SC(x, y) from SE(x, y) by EQ. 3. 
[55] The region signal intensity VC{x, y) can also be determined from EQ. 4, if 

the gain of the imaging system is known. That is: 

VC(x,y) = IC ° (x,y) EQ.17 
y GxSE(x,y) 

[56] If only the relative region signal intensity {VCr(x, y)) is required, it can be 

calculated without knowledge of the system gain. That is: 

ic°(x, y y 

VC R (x,y)= nu^l) EQ - 18 

[57] Figure 3 shows calibration standard 10 which has a top-left configuration. 

There are many other configurations that can be used to calculate shading errors and 
the region signal intensities. For example, Figure 8 shows a calibration standard 60 
with a bottom-right configuration. As with the top-left calibration standard 10 described 
in relation to Figure 3, this bottom-right calibration standard 60 has four detachable 
parts with N x M uniformly arranged calibration regions 5. The bottom-right corner 24 
consists of one calibration region 5 at the center; the bottom-row part 22 consists of N-1 
calibration regions 5; the right-column part 23 consists of M-1 calibration regions 5, and 
the part 21 consists of (N-1) x (M-1) calibration regions 5. The calibration regions 5 are 

referred to as R B/ W). 

[58] Figure 9 shows a horizontally shifted-and-rotated arrangement 61 of the 

bottom-right calibration standard 60. The shift-and-rotation is done by shifting the entire 
calibration standard 60 (in Figure 8) to right by one column and then moving the part 23 
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and 24 to the left of part 21 and 22. The re-arranged calibration regions 5 are referred 
to as ¥ HR (x,y), and: 

Content(R™ (x,yj) = Content(R** (Mod(x - l,N),y)) EQ. 1 9 

forx=1,...,N and y=1,...,M 

[59] Figure 1 0 shows a vertically shifted-and-rotated arrangement 62 of the 

bottom-right calibration standard 60. The shift-and-rotation is done by shifting the entire 
calibration standard 60 (in Figure 8) down by one row and then moving the part 22 and 
24 to the top of part 21 and 23. The re-arranged calibration regions 5 are referred to 
asP KZ W), and: 

Content(R vo (x,yj) = Content(R 5R (xMod{y - 1,M))) EQ. 20 

forx=1,...,N and y=1,...,M 

[60] T BR (x,y) , Y HR (x,y) and? VD (x,y) are the pixel sets corresponding to 

regions R BR (x,y), R HR (x,y) andR VD (x,y) respectively. Similar to the top-left 

configuration, it is required that the re-arranged calibration regions of the horizontally 
and vertically shifted-and-rotated arrangements 61 and 62 are aligned to the regions in 
the initial calibration standard 60. That is: 

V BR (x,y) = -p HR (x,y) = ¥ VD (x,y) RQ. 2 

forx=1,...,N and y=1 M 

[61] Similar to the deduction of EQ. 11 , the equation for the shading error 

SE L (x, y) of the bottom-right calibration standard 60 is: 

ICS M (x,y) = 2 x SE L (x,y) - SE L {Mod{x - l,N),y) - SE L (x,Mod(y - 1,M)) EQ. 21 

forx=1,...,N and y=1,...,M 

Where: ICS BR (x,y) = ICf (x,y) + IC[ D (x,y)-IC B L R (Mod(x-l,N),y)-IC B L R (x,Mod(y-l,M)) 

icl R (x, y )\s the log transformed image intensity of the region 
R BR (x,y) shown in Figure 8; 
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ic™(x,y)\s the log transformed image intensity of the region R HR (x,y) 
shown in Figure 9; and 

ic v L D (x,y) is the log transformed image intensity of the region R VD (x,y) 
shown in Figure 10. 

[62] Similar to the deduction of EQ. 12 (the matrix form of EQ. 1 1), the matrix 

form of EQ. 21 can be easily deduced. EQ. 21 has the same properties as those of EQ. 
1 1 , and the procedure for solving EQ. 21 is exactly the same as that for solving EQ. 11, 
except the coefficients of the matrix corresponding to EQ. 21 is different than those in 
EQ. 12. 

[63] In the previous discussion, the solution of shading error SE(x, y) from a 

top-left or a bottom-right calibration standard is derived from the image intensities of the 
initial calibration standard (10 or 60), a horizontally shifted-and-rotated arrangement (11 
or 61) and the vertically shifted-and-rotated arrangement (12 or 62). There are other 
configurations and arrangements, which can be used to derive the shading error SE(x, 
y) for an imaging system. 

[64] Reference is next made to Figure 11. In the previous discussion, the 

calibration standards (10 and 60) presented are rectangular. Calibration standards may 
be in other forms as well. For example, calibration standard 310 is a parallelogram. 
Calibration standard 310 has four parts 301, 302, 303 and 304. Each part has a 
number of calibration regions 305. In fact, a calibration standard can be of any form as 
long as the calibration regions of shifted-and-rotated arrangements can be aligned with 
the regions of the initial calibration standard. 

[65] The corner part (e.g., part 4 in Figure 3) of calibration standard 10 (Figure 

3) consists of only one calibration region 5. Referring to Figure 12, a calibration 
standard may alternatively have more than one region in the corner part. Figure 12 
shows a calibration standard 30 with a top-left configuration and with 2 calibration 
regions in the corner part 34. The top-row part 32 consists of N-2 calibration regions; 
the left-column part 33 consists of 2 x (M-1) calibration regions, and bottom-right part 31 
consists of (N-2) x (M-1) calibration regions. The horizontally shifted-and-rotated 
arrangement corresponding to calibration standard 30 will require shifting TWO 
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columns. Similar to the deduction of EQ. 1 1 , we can also derive an equation to solve 
the shading error SE{x, y) from the image intensities of the initial calibration standard, a 
horizontally shifted-and-rotated arrangement, and a vertically shifted-and-rotated 
arrangement. However, unlike EQ. 1 1 , solving this equation requires specifying TWO 
values of SE(x, y), e.g., SE(1, 1) and S£(2, 1) are given. Therefore, more than one 
initial value of the shading error SE(x, y) needs to be defined if the corner parts of a 
calibration standard consists of more than one calibration region. In order to obtain the 
correct shading correction coefficients for an imaging system using a calibration 
standard 30, the ratio of the shading errors for the two specified regions must be known 
a priori. 

[66] In certain applications, fabrication of a physically detachable calibration 

standard may not be convenient or possible. In this case, the shift-and-rotation can be 
achieved by performing two shifts of a calibration standard. Figure 13 shows calibration 
standard 70 with N x M calibration regions RS(x, y) with an imaginary top-left 
configuration. The dotted lines 81 and 82 in Figure 13 represent the imaginary divisions 
to the standard 70 into a top-left configuration, and imaginary parts 71, 72, 73 and 74 
correspond to part 1 , 2, 3 and 4 in Figure 3. It is to be noted that the use of the word 
"imaginary" here indicates only that calibration standard 70 cannot be conveniently 
detached to separate parts 71, 72, 73 and 74. 

[67] Figure 14 shows the process to obtain the image of a horizontally shifted- 

and-rotated arrangement 42 of calibration standard 70 (Figure 13). 
[68] First, calibration standard 70 is shifted to the left by 1 column so that line 

81 is aligned with the left edge of the image of calibration standard 70 in its original 
position. Image 40 is acquired with calibration standard 70 in this position. The shaded 
area of image 40 contains regions in the imaginary parts 71 and 72 of calibration 
standard 70. 

[69] Second, calibration standard 70 is shifted to the right by N-1 column 

(relative to its original position) so that line 81 is aligned with the right edge of the image 
of calibration standard 70 in its original position. Image 41 is acquired with calibration 
standard 70 in this position. The shaded area of image 41 contains regions in the 
imaginary parts 73 and 74 of calibration standard 70. 
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[70] Finally, images 40 and 41 are combined to form horizontally shifted-and- 

rotated arrangement 42. Areas 71 and 72 in horizontally shifted-and-rotated 
arrangement 42 are copied from the corresponding areas in image 40, and areas 73 
and 74 in horizontally shifted-and-rotated arrangement 42 are copied from the 
corresponding areas in image 41. The combined horizontally shifted-and-rotated 
arrangement 42 is essentially the same as if calibration standard 70 had been detached 
along line 81 and then re-assembled in a horizontally shifted-and-rotated arrangement 
and then imaged. 

[71] Figure 1 5 shows the process to obtain the image of a vertically shifted- 

and-rotated arrangement 52 of calibration standard 70. 

[72] First, calibration standard 70 is shifted up by 1 row so that line 82 is 

aligned with the top edge of the image of calibration standard 70 in its original position. 
Image 50 is acquired with calibration standard 70 in this position. The shaded area of 
image 50 contains regions in the imaginary parts 71 and 73 of calibration standard 70. 
[73] Secondly, calibration standard 70 down by M-1 rows (relative to its original 

position) so that line 82 is aligned with the bottom edge of the image of calibration 
standard 70 in its original position. Image 51 is then acquired with calibration standard 
in this position. The shaded area of image 51 contains regions in the imaginary parts 
72 and 74 of calibration standard 70. 

[74] Finally, the image 50 and 51 are combined to form vertically shifted-and- 

rotated arrangement 52. Areas 71 and 73 in vertically shifted-and-rotated arrangement 
52 are copied from the corresponding areas in image 50 and areas 72 and 74 in 
vertically shifted-and-rotated arrangement 52 are copied from the corresponding areas 
in image 51. The combined shifted-and-rotated arrangement 52 is essentially the same 
as if calibration standard 70 had been detached along line 82 and then re-assembled in 
a shifted-and-rotated calibration array and then imaged. 

[75] The shift operations in Figure 14 and Figure 15 can be accomplished by 

using either manual or computer-controlled XY positioning devices. It is preferable to 
have an accurate and reliable positioning device. 

[76] It is desirable that the regions RS(x, y) of horizontally and vertically 

shifted-and-rotated arrangements 42 and 52 align precisely with the regions RS(x, y) of 
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the initial calibration standard. Due to various mechanical errors, this precise alignment 
may not be possible, especially in a microscopy application. Using a calibration 
standard with a top-left configurations as an example, a mis-alignment is indicated by: 
P(x, y) * P H (x, y) or P(x, y) * P v (x, V) for some or all (x, y). For example, P(1 , 1) = 
{pxl(1, 1), pxl(2, 1), pxl(3, 1)} and P H (1, 1) = {pxl(2, 1), pxl(3, 1), pxl(4, 1)}. If P(x, y) is 
used as reference, i.e., the calibration regions of the initial setup of the calibration 
standard, the mis-alignment of P H (x, y) and/or P v (x> y) with P(x, y) introduces 
measurement error in IC H (x 9 y) and/or IC v (x,y). If the mis-alignment is significant, 
local correction functions (one for each calibration region) can be established and used 
to correct the measurement error in IC H (x 9 y) and/or IC v (x,y) due to the mis- 
alignment. Given image intensity /C(x, y) of a region with corresponding pixel set P'(x, 
y), this local correction function estimates the imaging intensity /C'(x, y) of the same 
region at pixel set P(x, y). This location correction function can be established by 
shifting and imaging the initial setup of a calibration standard around its default location. 
The image intensities of the calibration regions at different locations can then be used to 
define the local correction function parameters, using standard mathematical 
procedures. The established correction function should then be incorporated in the 
procedure 105 as shown in Figure 7 to correct the measurement error if mis-alignment 
happens. 

[77] The measured image intensities, such as IC°(x 9 y) , IC H (x 9 y) 

andIC v (x,y) in a top-left calibration standard (i.e. calibration standard 10 (Figure 3)), 

are not error free. In additional to the mechanical positioning error mentioned above, 
the other sources of error are electronic, quantum noises and etc. These measurement 
errors may cause errors in calculation of shading errors, shading correction coefficients 
and region signal intensities. 

[78] For a top-left calibration standard, IC°(x,y) , IC H (x 9 y) andicT(jt, j>) are 

the true image intensities (without measurement errors). The actual measurements can 
be expressed as: 

IC\x 9 y)^lc\x 9 y)x{\ + e%x,y))^e° A {x 9 y) 
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IC H (x, y) = IC H (x,y) x (1 + e H P (x, y )) + <f (x, y) EQ. 22 

/C v (x,y) = ^C V (^y)x(l + ^U,y)) + ^(^y) 

for x=1,...,N and y=1,...,M 

where: j>) ,e$ (x,y)ande v P (x,y) are the proportional measurement errors, 

and 

e A (x,y),e A (x,y)ande A (x,y) are the additive measurement errors. 
[79] The corresponding solutions of shading error can be expressed as: 

SE(x,y) = SE(x,y)x(\ + e SE (x,y)) EQ. 23 

forx=1,...,N and y=1 J ...,M 

where: SE{x,y) is the true shading error without error, and 

e SE (x,y) is the calculation error due to the measurement error 

e° P (x, y) , e H P (x, y) , e v P (x, y) , e° A (x, y) , (x ? y) and (x, 3;) . 
[80] The calculation errors (x,y) is summarized as a proportional error 

because of the relative nature of the shading error. Errore 5i? (x, j/) depends on the 
measurement errors and size of the calibration standard (i.e., N and M). The larger the 
measurement errors, the larger the calculation error e SE (x,y) , and the larger the size of 

the calibration array, the larger the calculation error e SE (x,y) . 

[81] Similarly, we can define the calculation errors for other calibration 

standards (bottom right, etc.). 

[82] The magnitude of errors can be reduced by averaging multiple 

acquisitions and calculations from the same calibration standard, by averaging the 
results of different calibration standards, or by averaging the combination of both. 
[83] As an example of averaging multiple calculation from the same calibration 

standard, SE (l) (x,y) an dSE (2) (x,y) ma y be defined as the results of two acquisitions and 

calculations of shading errors of an imaging system from a calibration standard with a 
top-left configuration. The averaged shading error can be defined as: 
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SEV(x 9 y) + SE«\x 9y ) 
SE(x, y) = EQ. 24 

[84] As an example of averaging different calibration standards, SE TL (x 9 y) may 

be defined as the results from a top-left calibration standard and SE BR (x 9 y) may be 

defined as the results from a bottom-right calibration standard. A simple averaged 
shading error can be defined as: 

SE Tl (x 9 y) + SE BR (x 9 y) 
SE(x 9 y) = Ll£^ \^IL EQ. 25 

[85] The statistics (such as mean and variance) of calculation error e SE (x 9 y) is 

position-dependent, even if such statistics of measurement error is position-independent 
(e.g., all e$(x 9 y) >ep (x 9 y) ,e v P (x 9 y) ,e° A (x 9 y) ,e A (x 9 y)ax\de v A (x 9 y)axe Gaussian noise with 
ji = 0 and u = 0.1). By using standard statistical methods and knowledge of 
measurement error, the error mean ^ese(x, y) and error standard deviation oese(x, y) of 
e SE (x 9 y) can be estimated for given measurement errors (e.g., Gaussian noise with ji = 
0 and <j= 1 for the additive measurement errors and Gaussian noise with /n = 0 and cr = 
0.05 for the proportional measurement errors). Using the error statistics ve$e(x, y) and 
0ese(x, y), weighting factors can be derived. Instead of using a simple averaging, a 
weighted averaging can be performed. 

c*v ^ ^ TL (^y)><SE TL (x 9 yHw BR (x 9 y)xSE BR (x 9 y) 

SE(x 9 y) = Tj -rz EQ. 26 

w TL (x 9 y) + w BR (x 9 y) 

where: w TL (x 9 y) and w BR (x 9 y) are the weighting factors for the top-left 

configuration and bottom-right configuration respectively. 
[86] The weighted average will reduce the overall error and will also reduce the 

variance of error statistics. The weighting factor can exist in many forms. For example: 



wTL (^y) = —nr^ r and wBR (^y)^~wr? r EQ. 27 



[87] Calculation error reduction of shading error will reduce the calculation 

error of shading correction coefficients. EQ. 23 to EQ. 27 describes error analysis and 



-23- 



averaging methods for shading error. Similarly, we could also directly analyze and 
average shading correction coefficients as well. 

[88] In calibration standard N and M were defined to be greater than 1 . If one 

of N or M is equal to one, the system becomes a 1 D (1 -dimensional) case. For 
example, the image sensor array consists of one row or one column of sensor 
elements, or only one row of regions is of interest. Without loss of generality, the 
following discussions will be based on a horizontal 1D array. A vertical 1D array will be 
simply rotated 90 degree. The naming convention for the calibration standards and 
shifted-and-rotated arrangements is similar to the 2D case. 

[89] Figure 16 illustrates a one dimensional calibration standard 100 with a left 

configuration. The calibration standard 100 contains N uniformly spaced regions. In the 
1D case, N has to be greater than one. The left part 102 consists of one region 5, and 
right part consists of N-1 regions 5. Figure 17 shows the spatial dimensions of standard 

100. "a" is the interval of regions 5. The regions 5 are referred to as Rl(x). 

[90] Figure 1 8 illustrates a shifted-and-rotated arrangement 1 10 of calibration 

standard 1 00. The shift-and-rotation is done by shifting the entire calibration standard 
100 (in Figure 16) to right by one column and then moving part 102 to the right of part 

101. The re-arranged calibration regions 5 are referred to asRl^OO , and: 

Content(Rl s * (jc)) = Content(Rl(MoJ(x + \,N))) for x= 1 , . . . , N EQ. 28 

[91] As with the 2D case, calibration regions of the shifted-and-rotated 

arrangement 110 are required to be in alignment with the regions 5 of the initial 
calibration standard 100. That is: 

Pl SR (x) = Pl(*) forx=1,...,N RQ.3 

where: pi(x) is the pixel set corresponding to region R1 (x) , and 

Pl ss O) is the pixel set corresponding to region R1 S *0) . 

Similar to the 2D case, the following relationship can be developed: 

ICl°(x) = GlxSEl(x)xVCl(x) forx=1,...,N EQ. 29 

where: lCl°(x)\s the image intensity of the region ri( X; y ) shown in Figure 16. 
SE1(x) is the shading error of the region Ri(x,j/) - 
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VC1{x) is the signal intensity of the region Rl(x, y) . 
G1 is the imaging system gain. 

And: 

ICl?(x) = Gl L +SEl L (x) + VCl L (x) forx=1,...,N EQ. 30 

IClf (jc) = G1 L + SE1 L (x) + VC1 L (Mod(x + 1,N)) for x= 1 N EQ. 31 

where: ici°{x)= In (ici°(x)) 

ici s L R {x)= In {ici SR {xj), and ici SR (x)\s the image intensity of region 5 

Rl Sfl (x, y) shown in Figure 18. 
SE1 L (x) = In (SE1 L (x)) 
VC1 L {x) = In (VC1(x)) 
G1 L =\n (G1) 

And: 

ICSl(x) = SEl L (x) - SEl L (Mod(x + 1, N)) for x=1 N EQ. 32 

where: /cs7(jc) = icif (x) - /Cif (Mod(x +i,N)) 
[92] As with the 2D case, SE1 L (x) can be resolved from EQ. 32 if one of the 

values is specified. 

[93] Similar to the 2D case, SE1(x) is defined as the solution (intermediate 

shading error) of EQ. 32 when SE1 L (N) = 1, and the shading errors, defined in EQ. 33, 
are uniquely defined. 



SEl(x) = —-if 



£ forx=1,...,N EQ.33 



e sim 



[94] Once the shading error is calculated, the shading correction coefficients 

SC1{x), signal intensity VC1(x), and relative signal intensityVC^x) can be calculated. 

SCl{x) = —1— for x=1 ,. . . ,N EQ. 34 

SEl(x) 

VCl(x)= IC1 ° (X) forx=1,...,N EQ.35 

GlxSEl(x) 
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ICl°(x)/ 

VCl R (x)= N /SE1( ) X) forx=1 N EQ.36 

1 V/C/W 

N f x /sEm 

[95] As with the 2D case (Figure 14), the shading errors can also calculated 

from a right configuration as well. Figure 16 and Figure 18 show that the right 
configuration is symmetrical to the left configuration. That is, the initial setup of a right 
calibration standard is the shifted-and-rotated arrangement 110 of a left calibration 
standard and vice versa. Therefore, there is no difference between a left configuration 
and a right configuration, except in the order they are given to an imaging system. 
[96] As with the 2D case (Figure 14), the shift-and-rotation can be achieved by 

two shifts in the case where the detachable configuration is not convenient or 
impossible to fabricate. 

[97] We have presented both 1 D and 2D cases. The same shift-and-rotation 

principle can be extended to the 3D (3-dimensional) domain as well. In the 3D domain, 
the field of view is 3 dimensional, and a calibration region becomes a calibration cube. 
The extension of a 2D calibration standard with a top-left configuration is a 3D 
calibration standard with a top-left-front configuration as shown in Figure 19. The 
calibration standard 400 has N x M x L calibration cubes and has 8 detachable parts. 
Part 401 consists of (N-1) x (M-1) x (L-1) calibration cubes; part 402 has (N-1) x (L-1) 
calibration cubes; part 403, which is not visible (blocked by parts 406 and 408) in this 
figure, consists of (M-1) x (L-1) calibration cubes; part 404 consists of 1 calibration 
cube; part 405 consists of (N-1) x (M-1) calibration cubes; part 406 consists of (L-1) 
calibration cubes; part 407 consists of (N-1) calibration cubes; and part 408 consists of 
(M-1) calibration cubes. The configuration is named as top-left-front because of the 
position of part 404 is at top-left-front. A horizontally shifted-and-rotated arrangement of 
the calibration standard 400 will be done by shifting the entire standard 400 to the left by 
one column and then moving parts 403, 404, 406 and 408 to the right side of the parts 
401, 407, 402 and 405 respectively. Similarly, the vertically shifted-and-rotated 
arrangement and the front-to-back shifted-and-rotated arrangement will be form in a 
similar manner. 
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[98] The shading errors SE3{x, y, z) of the calibration cubes can be resolved 

from the 3D image intensities of the initial setup, the horizontally shifted-and-rotated 
arrangement, vertically shifted-and-rotated arrangement and the front-to-back shifted- 
and-rotated arrangement. All the mathematical deductions are similar to the 2D case. 
[99] Reference is again made to Figure 3. The method of the present 

invention allows the log-transformed shading error {SE L in EQ. 11 or EQ. 12) of each 
region 5 to be represented as a factor of the log-transformed shading errors of other 
regions 5. An intermediate shading error for each region may be calculated using 
standard matrix techniques, if the shading error for one region is specified (assuming 
that corner part 4 has a single region so that matrix A in EQ. 12 has a rank of N x M -1). 
The shading error may then be calculated as the normalized value of the intermediate 
shading error (EQ. 13 or EQ. 16). The effect of the selection in the intermediate 
shading errors is eliminated by the normalization. 

[100] Reference is made to EQ. 13, EQ. 16 and EQ. 33. Due to relative nature 

of shading correction, the normalization to derive a unique shading error is not 
completely necessary. In other words, the intermediate shading errors could be used to 
derive the shading correction coefficients. In this case, there may be an embedded 
system gain to the subsequently corrected images. 

[101] Referring to Figures 5 and 6, calibration standard 10 is re-arranged into 

horizontally shifted-and-rotated arrangements 11 and 12. Any other re-arrangements of 
calibration standard 10 which (i) allow the shading error of each region 5 in calibration 
standard 10 to be defined in terms of the shading errors of other regions and (ii) allow 
the shading error of each region 5 to be solved by specifying an shading error for one 
(or more) regions, may be used in the method of the present invention. A person skilled 
in the art will be capable of developing equations corresponding to EQ. 12 which allow 
the relative shading error for each region 5 to be solved. 

[102] Several embodiments have been described in details as examples of the 

present invention. These embodiments may be modified and varied without departing 
from the spirit and scope of the invention, which is limited only by the appended claims. 
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